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Abstract 

Power transformer outages have a considerable economic 
impact on the operation of an electrical network. In order to 
draw maximum power from transformers and, at the same 
time, avoid thermal mishaps, it is essential to carefully study 
its thermal behavior. Furthermore, an accurate computation 
of the hottest spot temperature (HST) helps in a realistic 
estimation of the reliability and remaining life of the 
transformer winding insulation. This paper presents steady 
state temperature distribution of a power transformer layer- 
type winding using conjugated heat transfer analysis, 
therefore energy and Navier-Stokes equations are solved 
using finite difference method. Meanwhile, the effects of 
load conditions and type of oil on HST are investigated 
using the model. Oil in the transformer is assumed nearly 
incompressible and oil parameters such as thermal 
conductivity, special heat, viscosity, and density vary with 
temperature. Comparing the results with those obtained 
from finite integral transform checks the validity and 
accuracy of the proposed method. 
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Introduction 

In a power transformer a part of the electrical energy 
is converted into the heat. Although this part is quite 
small comparing to total electric power transferred 
through a transformer, it causes significant 
temperature rise, which represents the limiting criteria 
for possible power transfer through a transformer. 
That is why the precise calculation of temperatures in 
critical points (top oil and the hottest solid insulation 
spot) is of practical interest. Thermal impact leads not 
only to long-term oil/paper-insulation degradation; it 
is also a limiting factor for the transformer operation 
[1], [2]. Therefore, the knowledge of the temperature, 
especially the hottest spot temperature, is of high 



interest. If the temperature rise goes beyond the 
permissible value, the load of transformer must be 
reduced or an auxiliary transformer is used in order to 
preserve the insulation from deterioration. For an oil- 
immersed transformer, the oil surrounds the 
transformer body. Oil is a nearly incompressible fluid 
and density changes due to temperature rise, therefore 
oil moves in the transformer. The heat transferred by 
convection is the most important method of heat 
transfer. 

Hottest spot temperature must not exceed the 
prescribed value in order to avoid insulation faults. A 
hottest spot temperature calculation is given in the 
International Standards [3]-[5]. The algorithm for 
calculating the hottest spot temperature of a directly 
loaded transformer using data obtained in a short 
circuit heating test is given in [6] . Heat transfer theory 
results from winding to oil are exposed in [7], [8]. 

In this paper, author has proposed a procedure for 
obtaining the temperature distribution in the power 
transformer and the effects of load conditions and type 
of oil are investigated using the model. For this reason 
energy and Navier-Stokes equations are solved using 
finite difference method. Therefore, a code has been 
provided under MATLAB software. The model can be 
used for temperature calculation on the arbitrary 
change of current and outside air temperature. In the 
paper, thermal model, energy, and Navier-Stokes 
equations are given in section 2. Results and 
discussion of the proposed work have been provided 
in section 3. 

Thermal Model 

Energy equation for Newtonian incompressible fluid 
(or nearly incompressible fluid) such as oil in two 
dimensions is [9] 
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d(C-T) d(C-T) 
p ■ (V — ^ + V — ^ -) = 
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dx dx By By 



Continuity Equation: 



(1) 



Where T is temperature, Vx is velocity in x direction, V y 
is velocity in y direction, k is thermal conductivity, p 
is density and Q is special heat of the oil in the 
transformer. In equation (1), oil properties vary with 
temperature [1], [10]-[12]. The temperature 
dependence of oil properties are given as follows: 



H = a x - exp( 



1 + 273 



) 



C p = a 3 + a 4 -T 



p = a 5 +a 6 T 
k = a 1 +a$-T 
P = a, 



(2) 

(3) 
(4) 
(5) 
(6) 



Where//: is viscosity and p is volumetric expansion 
coefficient. The nine constants for transformer oil and 
silicon oil have been listed in Table 1 [3, 4]. It is 
generally valid for all types of transformer oil that the 
variation of the oil viscosity with temperature is much 
higher than the variation of other oil properties [11], 
[12]. Thus, all oil physical properties except the 
viscosity can be replaced by a constant. However, the 
variations of all oil properties with temperature have 
been considered in this paper. 

TABLE 1 OIL CONSTANTS 



Oil constant 


Transformer oil 


Silicon oil 


ai 


0.0000013573 


0.00012127 


az 


2797.3 


1782.3 


as 


1960 


1424 


at 


4.005 


2.513 


05 


887 


989 


06 


-0.659 


-0.870 


07 


0.124 


0.138 


as 


-0.0001525 


-0.00009621 


09 


0.00086 


0.00095 



In (1) velocity is unknown, and then we must solve 
Navier-Stokes equations. Navier-Stokes equations in 
two dimensions for Newtonian incompressible fluid 
are [9] 



dV. d V, 



d x 



+ ■ 



= 



d y 

Momentum Equations: 



8V 8V dP 

p.(V x — + v ,—) = - — 

dx y dy dx 

d dV d dV dV v 

+ (—(2- M -^) + —( M -^ + M -^)) 

dx dx dy dy dx 



(7) 



(8) 



dV dV dP 
p.(V x ^ + V v ^) = - — + 
dx ' dy dy 

d dV dV v d dV y 

(—(M- — + M-—) + — (2-M-—)) 

dx dy dx dy dy 

+ p . g ./3.(T-T ref ) 



(9) 



Where P is pressure and g is gravitational 
acceleration .The heat conduction equation for core 
and windings is written [10]-[12] 



d 2 T . d 2 r n n 

d x 2 y d y 2 



(10) 



Where kx and k y are thermal conductivity in x and y 
directions respectively. The term Q is the volumetric 
heat source function and has been modified here to 
take care of variation of electrical resistance of copper 
with temperature. The heat source term Q can be of 
the form 



Q = Qo-\X + a e <T-T^y] 



(11) 



Where ac is the temperature coefficient of electrical 
resistance of copper wire. With this representation, the 
function Q becomes temperature dependent, 
distributed heat source. Thermal conductivities are 
unequal in different directions. Thermal conductivity 
has been treated as a vector quantity, having 
components in both radial and axial direction. 
Resultant thermal conductivity of the system is [10] 



K 



(12) 
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Where 



log^ 



1 r 2 1 ^ 

log — log — 
r, r 2 

+ - + ... + • 



log^ 



£ _ ' ' ( f cu + tin ) 

y t- -k +t -k. 

in cu cu m 

Term K represents resultant thermal conductivity of 
insulation and conductor system. In this paper, the 
bottom oil temperature rise over ambient temperature 
has been calculated as 



0„ 



R, +1 



(13) 



Where 0u is the bottom oil temperature rise over 
ambient temperature, 6p is the full load bottom oil 
temperature rise over ambient temperature obtained 
from an off-line test; Ri is the ratio of load loss at rated 
load to no-load loss. The variable L is the ratio of the 
specified load to rated load. 



(14) 



ruled 



The exponent n depends upon the cooling state. The 
loading guide recommends the use of n=0.8 for natural 
convection and n=0. 9-1.0 for forced cooling. 

Results and Discussion 

A code has been provided for finite difference solution 
of energy and Navier-Stokes equations using 
MATLAB 7.1. Figure 1 shows cross section of a 50KVA, 
20KV/400V power transformer in two dimensions. 
Dimensions and specifications of the power 
transformer have been summarized in Table 2. Table 3 
shows the power transformer losses. 

We assume that inlet oil temperature is the 
temperature base. Therefore, we have 



T -T 
AT = ^xlOO 



T 



(15) 




FIGURE 1 THREE PHASE POWER TRANSFORMER IN 
TWO DIMENSIONS 

TABLE 2 DIMENSIONS AND SPECIFICATIONS OF THE 
PROPOSED POWER TRANSFORMER 



Core diameter (cm) 


10 


Width of each window (cm) 


12 


Height of window (cm) 


28 


Thickness of LV winding (cm) 


1.2 


Thickness of HV winding (cm) 


2 


Rated power 


50KVA 


HV voltage 


20kV 


LV voltage 


400V 


HV current 


1.44 A 


LV current 


72A 


TABLE 3 POWER TRANSFORMER LOSSES 


Losses (w) 


Value 


Core 


158 


DC of LV windings 


384 


DC of HV windings 


534 


Eddy currents of LV windings 


64 


Eddy currents of HV windings 





Where Tb,on is the bottom oil temperature. Figure 2 
shows temperature distribution from LV winding 
(axis symmetry in height direction) in one per unit 
(p.u.) load with oil natural cooling (ON). It can be 
pointed out that for LV winding, maximum 
temperature location is around 80% of winding height 
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from the bottom and at about 50% of radial thickness 
of the layer. Temperature distribution from HV 
winding has been shown in figure 3. It can be 
observed that the maximum temperature occurs in the 
neighborhood of 55% of the axial and 50% of the radial 
thickness of the layer. 
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FIGURE 2 TEMPERATURE DISTRIBUTIONS IN HEIGHT 
DIRECTION OF LV WINDING WITH TRANSFORMER OIL 
NATURAL COOLING 
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FIGURE 3 TEMPERATURE DISTRIBUTIONS IN HEIGHT 
DIRECTION OF HV WINDING WITH TRANSFORMER 
OIL NATURAL COOLING 
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FIGURE 5 HST VERSUS LOAD AT Tamb=25(°C) WITH OIL 
NATURAL COOLING 



Figure 4 shows HST rise over bottom oil temperature 
versus load for two types of oils and figure 5 shows 
HST versus load at Tamb=25(°C). 
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FIGURE 4 HST RISE OVER BOTTOM OIL TEMPERATURE 
VERSUS LOAD AT T™i,=25(°C) WITH OIL NATURAL 
COOLING 

It can be pointed out from figures 4 and 5 that the HST 
value is low when the transformer oil is used for 
cooling, therefore the transformer oil is better than 
silicon oil in the cooling of the power transformer 
when load changes. Table 4 shows comparison of the 
proposed method with finite integral transform used 
in [10] for power transformer oil. HST locations at 
different loading are given in Table 5. 

TABLE 4 HST (°C) MAGNITUDES AT T™^25(°C) 



Load (p.u.) 


Proposed 


Analytical [10] 


0.6 


51.72 


53 


0.7 


58.35 


56 


0.8 


65.71 


66 


0.9 


73.78 


75 


1.0 


82.54 


80 


1.1 


91.99 


93 


1.2 


102.15 


101 


1.3 


113.01 


115 


1.4 


124.59 


125 


1.5 


136.91 


139 



Reference [10] has proposed an analytical method 
based on boundary value problem of heat conduction 
in power transformer winding using finite integral 
transform techniques. This technique uses Fourier 
transform and Hankel transform and finds 
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temperature based on polynomials with infinite order 
(sequential series) of space variables and time, and 
then approximates temperature with finite order. 
Analytical method in [10] is suitable for cylindrical 
winding with rectangle cross-section, but this 
limitation is not for finite difference method. In the 
other hand, if we assume that we can increase order of 
polynomial, solution time will rise and analytical 
method will not differ with numerical method. 



TABLE 5 HST LOCATIONS AT T Ilm6 =25(°C) 



Load (p.u.) 


Proposed (radial 
thickness%-height%) 


Analytical [10] (radial 
thickness%- height%) 


0.6 


50 / 54.17 


44/52 


0.7 


50 / 54.17 


45 / 53.5 


0.8 


50 / 54.58 


46/54 


0.9 


50 / 54.58 


48 / 54.5 


1.0 


50/55 


50/55 


1.1 


50/55 


50/55 


1.2 


50 / 55.42 


50/55 


1.3 


50 / 55.42 


50/55 


1.4 


50 / 55.42 


50 / 55.5 


1.5 


50 / 55.42 


50 / 55.5 



Conclusion 

A numerical study was conducted to investigate 
temperature distribution in layer-type power 
transformer for a variety of load conditions and type 
of oil. Heat transfer partial differential equations were 
solved numerically with finite difference method. The 
heat source function has been taken as temperature 
dependent and has been directly incorporated in the 
heat conduction equation for core and windings. The 
purely numerical approach followed in this paper 
seems to correspond reasonably well with results of 
analytical calculations [10]. In addition, it can be 
pointed out that transformer oil is better than silicon 
oil in cooling of the power transformer when load 
changes. 
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Nomenclature 

C P : specific heat of oil (J.kg 4 .K 4 ) 



g: gravitational acceleration (m.s 2 ) 

L: ratio of the specified load to rated load 

k: thermal conductivity of oil (W.m 4 .K 4 ) 

K: equivalent thermal conductivity (W.m^.K 1 ) 

kcu : thermal conductivity of copper (W-m^.K 1 ) 

km : thermal conductivity of insulation (W.m lK 4 ) 

kx : thermal conductivity in x direction (W.m 4 .K 4 ) 

ky : thermal conductivity in y direction (W.m 4 .K 4 ) 

n : temperature rise exponent due to bottom oil 

P: pressure (N.nv 2 ) 

Q : volumetric heat source function (W.nv 3 ) 

Qo: reference volumetric heat source (W.nv 3 ) 

n's : radius of insulation and conductor layers (m) 

Ri : loss ratio = load loss/no load loss 

T: temperature (K) 

Tref : reference temperature (K) 

Tb,oi : bottom oil temperature (K) 

tcu. thickness of conductor (m) 

tin : thickness of insulation (m) 

Vx. velocity in x direction (m.s 4 ) 

V y : velocity in y direction (m.s 1 ) 

Abbreviations 

HST: hottest spot temperature 

HV: high voltage 

LV: low voltage 

ON: oil natural cooling 

p.u.: per unit 

Greek symbols 

p : coefficient of volumetric expansion of oil (K 4 ) 

Bu : bottom oil temperature rise over ambient 
temperature (K) 

6fi: bottom oil temperature rise over ambient 
temperature at rated load (K) 

p : density of oil (kg.nv 3 ) 

po: reference density of oil (kg.nv 3 ) 

«r : temperature coefficient of electrical resistance (K 4 ) 

p. : viscosity of oil (kg.nv 1 . s 4 ) 
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